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Abstract 

Gravitational waves act like lenses for the light propagating through them. 
This phenomenon is described using the vector formalism employed for ordinary 
gravitational lenses, which was proved to be applicable also to a non-stationary 
spacetime, with the appropriate modifications. In order to have multiple imaging 
an approximate condition analogous to that for ordinary gravitational lenses must 
be satisfied. Certain astrophysical sources of gravitational waves satisfy this con- 
dition, while the gravitational wave background, on average, does not. Multiple 
imaging by gravitational waves is, in principle, possible, but the probability of 
observing such a phenomenon is extremely low. 
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1 Introduction 



In general relativity (and in all metric theories of gravity) a gravitational field deflects 
light rays propagating through it, and gravitational waves are no exception. It is well 
known that exact solutions of the Einstein equations representing plane gravitational 
waves exhibit focusing properties on timelike and null geodesies propagating through 
them PJ. Light propagation through realistic linearized gravitational waves outside the 
laboratory has also been studied [|]-||12|. Most authors restricted their attention to 
cosmological gravitational waves; a particular attention was payed to the frequency shift 
effect induced by cosmological gravitational waves on the microwave background pho- 
tons [TJ], H [T3|, |T5| , expecially after the COBE discovery of anisotropies in the cosmic 
microwave background. Other works focussed upon the effects of gravitational waves 
perturbing ordinary gravitational lenses [16|-[18|. In recent papers, the case of gravita- 



tional waves generated by astrophysical sources was studied fl9fl-f22|| and, in some cases, 
optimistic statements were made about the detectability of these waves via the effects 
induced by their interaction with the light coming from distant objects. 

The study of gravitational lensing by mass concentrations perturbing the background 
curvature of the universe has become a very active field of astronomy 111 extreme 
conditions multiple images of a celestial object are created, which are accompanied by 
high amplification events. In these situations spectacular phenomena are observed by 
the astronomers: multiple quasars, giant arcs, arclets and radio rings [ 23fl . It is therefore 
natural to ask whether gravitational waves can create multiple images of a celestial ob- 
ject and give rise to the associated high amplification events. This question was posed in 
Refs. 



24 - 27 



The possibility of multiple imaging by gravitational waves in general rel- 
ativity constitutes the subject of the present paper. Some of the ideas exposed here were 
anticipated in a brief abstract |27[] . While the idea of multiple imaging by gravitational 
waves was only qualitatively sketched in Ref. p4" |, in Ref. p3 the Raychaudhuri equa- 
tion was used to deduce that the amplification of a light beam induced by gravitational 
waves is of second order in the wave amplitudes, and therefore completely negligible 
in any situation of astrophysical interest. This conclusion is valid only when a single 
image of the light source is created, and cannot be extended to the case of multiple im- 
ages and the associated high amplification events because the scalar formalism is known 
to fail in these situations p8| . Nowadays, we have formalisms that are more suitable 
than the Raychaudhuri equation for the description of multiple images in gravitational 
lensing f23fl . There is, however, a problem with the available formalisms, i.e. they were 
introduced to describe lensing in a (conformally) static spacetime and, in principle, their 
validity is restricted to this class of spacetimes. Instead, the case considered in this paper 
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involves a highly non-(conformally) stationary spacetime. This problem was resolved 



in a previous paper [pfiH , where it was proved that the vector formalism |23| can still be 
applied to the case of lensing by gravitational waves, provided that a new formula for 
the deflection angle is calculated. This result may look intuitive, but the proof of its 
validity is non-trivial and requires a generalization of the Fermat principle to the case 
of non-stationary spacetimes, which was given only recently . 



Using the results of Ref. ||29| , we approach the problem employing the vector for- 
malism and describing the effect of the lens with a plane-to-plane map. The vanishing 
of the Jacobian of this geometric map is associated with the failure of the inverse map, 
and defines critical lines and caustics. An approximate condition for multiple imaging 
is derived, and it is shown to be analogous to the well known approximate condition for 
multiple imaging by ordinary gravitational lenses pi] . Order of magnitude estimates 
show that the condition for the creation of multiple images by gravitational waves is 
satisfied by certain astrophysical sources of waves considered in the literature. This 
happens thanks to the balance between the (large) source-lens distances and the (small) 
wave amplitudes, and because the lensing waves in the sky can have relatively large 
amplitudes (i.e. larger than those expected in the Solar System). On average, the 
gravitational wave background is not capable of multiple imaging. It is concluded that 
multiple imaging by gravitational waves is, in principle, possible, but the probability of 
observing such a phenomenon is very low. 

Although various aspects of lensing could be considered, in this paper we focus on 
the creation of multiple images because the associated high amplification events are the 
most dramatic effects, and the most relevant for the detection of lensing events and, 
indirectly, of gravitational waves. 

The plan of the paper is as follows: in Sec. 2 and Sec. 3 we introduce the necessary 
tools and the assumptions underlying this work. Emphasis is given to the computation 
of the deflection angle (Sec. 2) and to the validity of the thin lens approximation (Sec. 3). 
In Sec. 4, the modified vector formalism is applied to lensing gravitational waves and 
an approximate condition for the creation of multiple images is derived. In Sec. 5, 
gravitational waves and ordinary gravitational lenses are compared, while Sec. 6 explores 
the possibility that the condition for multiple imaging is satisfied by astrophysical sources 
of gravitational waves. Sec. 7 contains final remarks. 

We use units in which the Newton constant G and the speed of light c assume the 
value unity (but G and c will occasionally be restored). The metric signature is +2; 
Latin indices run from 1 to 3, and Greek indices run from to 3. A comma denotes 
partial differentiation. 
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2 Computation of the deflection angle 

Let us consider gravitational waves localized in a region of space between a light source 
and an observer. The spacetime metric is given, in an asymptotically Cartesian coordi- 
nate system {t, x, y, z}, by 

where rj^ =diag(— 1, 1, 1, 1) and \h^ u \ <C 1 in these coordinates. Let us consider a light 
ray whose unperturbed path is parallel to the z-axis. The photon describing this ray 
has four-momentum 

p" = PU + fc" = (1 + 8p°, Sp\ Sp 2 , 1 + 5p 3 ) , (2.2) 

where 5p^ are small deflections (of order h pv ) and the unperturbed photon four-momentum 
is p^ = (1,0,0,1). The computations are performed in the geometric optics approxi- 
mation, that holds if the wavelength X gw of the gravitational wave is much larger than 
the photon wavelength A em , and if 

Aem > Xgw 2 /Di , (2-3) 



where is the observer-lens distance. Equation ( |2.3j ) guarantees that the size of the 
interference fringes which eventually form at the observer's position is not comparable 
to the "geometrical shadow" of the lens [25]. In order to make computations to second 
order in the wave amplitudes, tensor indices are lowered and raised with and with 

g iw = V ^ _ h ^ + (h 2 ) = rf v - ifrj^h^ + 0{h 2 ) , (2.4) 

respectively. The symbols 0(h), 0(h 2 ) denote the orders of magnitude of h pu and of 
terms quadratic in h^ v , respectively. The equation of null geodesies gives 



d(5p^ 



+ (p p m + dp") (pfo + 5p°) = , (2.5) 



dX 

where A is an affine parameter along the null geodesies and 



K° = \ 9^ + hua, P - h pa ,„) . (2.6) 



One obtains 



d(5pf 1 ) ( 11, 

+ if" [ h u0 ,o + h u0 ,3 + h u3t0 + h v3t3 - h 03 , u - - h 00jV - - h 33 ,u ) + 0(h 2 ) = . 



(2.7) 
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We introduce the notation 



A x =p"A a ; B^p*—. (2.8) 

Using the expression 

T^a = VT [h v x,x - \ hxx,u) + 0(h 2 ) (2.9) 
and Eq. ( |2.7|) , one obtains 

W = ~ J° d\ (h\ x ~ \ hxx^ + 0(h 2 ) , (2.10) 

where the indices are now raised with jf" and the integral is computed along the photon 
path from the source to the observer. Since a localized pulse of gravitational waves 
is considered, the first term in the integrand gives a vanishing boundary termQ and it 
follows that 

W= \ 1° d\ (h 00 + 2h 03 + h 33 )'» + 0{h 2 ) . (2.11) 
2 Js 

Performing the integration along the unperturbed photon path instead of the actual path 
involves only an error of order h 2 and hence the first order result is given by 

Srfq = \ f s dz (h 00 + 2h 03 + h 33 y» (2.12) 

(notice that the deflection angle is the opposite of the angle used in the theory of 
ordinary gravitational lenses, in which deflection from a straight line is defined as having 
the opposite orientation). For pulses of gravitational waves this integral is extended to 
a region of space of order P (the characteristic period of the waves) where the pulse is 
localized so that, in order of magnitude, 5p ~ P ■ h/P = h. If the gravitational wave 
background is considered, it might appear that the computation of the integral along 
the whole photon path from the source to the observer gives 5p ~ Dh/P. However, this 
secular effect is absent, i.e. the deflection does not cumulate with the travelled distance, 
due to the transversality of gravitational waves and to the equality between the speed of 
light and that of the random inhomogeneities of the medium (the gravitational waves), 
as shown in Refs. @, |, § §, [TT], 0. 

1 The boundary term is usually negligible [ p2| . The case of "pulses with memory" was considered in 
Ref. M. 
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We also compute the second order correction 5p%\ to the deflection, which was not 
given in previous references and will be used in Sec. 4. 8p%\ is obtained by inserting 
Eq. ( |2.12| ) into the equation of null geodesies 



dX 



T^ + 2T^5p'( 1) +O(h 3 ) = 



(2.13) 



where 8p fl = Sp^ + 8p%^ , which gives 



(2) 



rO 
IS 



dz 



- {dX + W + doK + dsK - d^h 0a - d»h Ba ) 5p 



(i) 



+/''''' ( h a0 fi + h a Q t 3 + h a3 fi + h a - i3 — h 03a — —hoo, a — 2^ 33 > a 



+ 0(h 3 ), (2.14) 



with 8p°i) given by Eq. ( |2.12|) and where the integration is performed again along the 
unperturbed photon's path instead of the actual path, the difference contributing only 
by a third order term. Equation ( |2.14| ) allows the computation of the second order 
contribution to the divergence 0a(8p a ) (the first order contribution being zero j25| ): 

d A (5p A ) = / dz \h\ A + h\ aA + h\ M + h A aM 

J s 



-d A (5p a {1) ) (h\ a + h\ a + h A ufi + h A a , 3 - h 0a 



A , ,A 

— ^3(7 



+ h Aa ^A ( h a o,0 + ^a0,3 + /Iq-3,0 + h a 3,3 ~ ^03,a 
+ h Aa (h a QflA + ^aO,3A + h a 3flA + h a3 ^A — ^03,aA 

where A = 1, 2. 



h 1 h 

2 "'00,0: — 2 33 '° 

-^-oo.qA — -^ 33 > aA 



(2.15) 



3 The thin lens approximation 

In the next section we will apply the vector formalism for gravitational lenses, with the 
appropriate modifications, to the case of lensing by gravitational waves. One of the 
assumptions underlying the vector formalism is the validity of the thin lens approxima- 
tion, i.e. that the size of the lens be negligible in comparison to the length of the path 
travelled by the photon. This approximation requires special care in the case of lensing 
gravitational waves, and is dealt with in what follows. 
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To be specific, we restrict to the case of gravitational waves generated by astro- 
physical sources and propagating with wavefronts that are approximately spherical at 
a sufficiently large distance from the source. In this situation, a photon coming from 
infinity and entering the gravitational wave will never leave it, hence the thin lens ap- 
proximation might seem to be inapplicable. However, the photon is appreciably deflected 
only during a small fraction of the time spent inside the gravitational wave, at the min- 
imum distance from the source of gravitational waves (where the wave amplitude is 
maximum). The region of space in which the photon is appreciably deflected has a size 
that is much smaller than the total length of the path travelled by the photon. In fact, 
let h ~ hsRs/r be the order of magnitude of the gravitational wave amplitude at a 
distance r from the source of waves, where Rs and hs are, respectively, the amplitude 
at the source and the Schwarzschild radius of the object generating gravitational waves. 
One has hs — eM/R, where e, M and R are the efficiency of the generation process, 
the mass and radius of the object, respectively. For neutron stars, M/R ~ 1/20, and 
e ~ 0.1 is a very generous estimate of the efficiency (see Sec. 6). These numbers give 
a deflection angle 5 ~ h ~ 5 ■ 10~ 3 (1 km/r). The minimum detectable separation 
between multiple images is of the order 5q ~ 10~ 4 arcseconds, marginally accessible by 
VLBI. The requirement 5 > 5q gives r ~ hsRs/S < 3 ■ 10 7 km. This distance is very 
small in comparison with the path travelled by the photon, usually of the order of many 
megaparsecs. Allowing a separation between the images of a light source smaller than 
5q ~ 10~ 4 arcseconds would still give high amplification events although it would be 
impossible to resolve the multiple images (this occurs, e.g. in microlensing by stars or 
compact objects ||23|| ). Even in this case the size of the region where the deflection takes 
place is negligible in comparison to the size of the photon path. 

A second argument supporting the validity of the thin lens approximation is the 
following: Due to the transversality of gravitational waves, a photon is not deflected by 
a gravitational wave propagating parallely to the photon path. To be more specific, let 
k be the (3-dimensional) wave vector of a gravitational wave of finite (short) duration. 
If the component k 3 is zero (i.e. the gravitational wave propagates orthogonally to the 
unperturbed photon path), it is clear that the gravitational wave is confined to a plane 
and that the thin lens approximation holds. One may conjecture that this approximation 
fails to be satisfied as \k 3 \ increases, but this would be incorrect. We give an explicit 
proof for a single monochromatic gravitational wave. In the extreme case \k\ = \k 3 \ 
(gravitational wave antiparallel[| to the photon path) the light rays are not deflected, 

2 If the gravitational wave propagates parallely to the light rays, the light never reaches the gravita- 
tional wave, due to the equality between the speed of gravitational waves c gw and the speed of light. 
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due to the transversality of gravitational waves. This fact can be deduced by Refs. |2], [IT] 
or by Eq. ( p. 12 ). In fact, the use of the transverse-traceless (TT) gauge in Eq. ( 2.12 ) 



gives /loo^ = ^03 T ' ) = ^33 = 0) which implies 5p A = 0. The vanishing of the deflection 
angle (a gauge-invariant quantity) in the TT gauge implies its vanishing in any gauge. 

If the direction of propagation of the gravitational wave and the (unperturbed) pho- 
ton path are not parallel nor orthogonal, the light rays can be approximated by zig-zag 
paths with (3-dimensional) unit tangent vector n, and the deflection depends on the 
quantity (h OQ + 2hoin l + hijn % n^) , which replaces the argument of the integral in Eq. (|2.12|) 



The vector n is decomposed into its components parallel and perpendicular to k: 

ft = ny + n± , (3.1) 
where ny = k(n- k/k 2 ^j and n± ■ k = 0. By introducing the quantities 

h \\=h 0i n\ , h 0± = h 0i n l ± , (3.2) 

h\\ || = hij nj| nl , h\\ j_ = hij riy n J ± , h±± = hij n l ± n J ± , (3-3) 

one can write 

h 00 + 2h 0i ri + hij n l n j = h 00 + 2 f/i y + h ± + h±\\j + h± ± + h\\ y . (3.4) 

The monochromatic gravitational wave is best described in the TT gauge in which 
ho/j, = h\\ y = h\\ j_ = 0. In this gauge 

/i o + 2 (/io || + h ± + h± y) +hx± + h\\\\ =h±± . (3.5) 

Therefore, the "component" of the gravitational wave propagating parallely to the (un- 
perturbed) photon path does not affect the light rays, and the thin lens approximation 
is valid. 

4 The vector formalism and the occurrence of mul- 
tiple images 

The vector formalism familiar from the description of ordinary gravitational lenses can 
be employed also for the case in which the lens is a gravitational wave, provided that the 

This situation may change in alternative theories of gravity in which c gw ^ c |p3| . 



7 



deflection angle ( 2.12 ) is used. The proof of this statement was presented in a previous 



paper |29| , and relies on a rigorous formulation of the Fermat principle valid for arbitrary 
(non-stationary) spacetimes which became available only recently |5U] . 



We set the geometry of the problem as customary in gravitational lens theory: we 
consider only waves which are localized near a plane z = z w to satisfy the thin lens 
approximation; the deflection takes place in this plane {lens or image plane in the usual 
language) [p3f| . This description is clearly adequate for the case of a gravitational wave 
propagating in a direction orthogonal to the z-axis. Naively, one might expect that 
the description becomes less and less correct as the direction of propagation of the 
gravitational wave becomes closer and closer to the z-axis, but this is not the case, as 
discussed in the previous section. 

Let x = (x,y) be the apparent source position in the lens plane, and s = (s x ,s y ) 
be the true source position (i.e. its position were the lensing wave absent). The plane 
orthogonal to the 2-axis and passing through the source is referred to as the source 
plane. The true source position in the lens plane is related to that in the source plane 
by identifying the two planes). The deflection is described by a two-dimensional vector 
field 5p A (x) (A = 1,2) in the lens plane. The action of the lens is described by a 
plane-to-plane mapping x A i — > s A , where s is given by the lens equation 

s A = x A + DlD^ (4i) 

and where Dl, Dl$ and Ds are the observer-lens, lens-source, and observer-source 
distances, respectively. As customary in gravitational lens theory, we can fit cosmology 
into the model by assuming the .D's to be angular diameter distances in a Friedmann- 



Lemaitre-Robertson- Walker (FLRW) universe | 23| . However, for the sake of simplicity, 



we assume that the background is flat and that the D } s denote Euclidean distances, with 



Ds = Dl + Dls- As described in Ref. fl29 |, the use of the Fermat principle generalized 



to non-stationary spacetimes [|30] allows one to derive the lens equation ( |4.1| ), with the 
deflection angle 5p A given by Eq. (|2.12|) . 

The map described by Eq. ([4.1|) has the Jacobian matrix 



j(s\_ (d^\_ (l + Dd x (Sp*) Dd y (5p*) \ 

\x) \dx B ) \ Dd x (6 p y) i + Dd y (6py))> l4 - /j 

where D = D L D LS /D S - The inverse matrix A = J -1 represents the amplification 
tensor, while its determinant A = Det(J) -1 is the (scalar) amplification. Since the 
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surface brightness is conserved during lensing whenever geometric optics holds []32 



and 



^ area of an infinitesimal region in the lens plane ^. 
area of the corresponding region in the source plane ' 



A has also the meaning of the ratio of light intensities with and without the lens [23 



A small circular source will be imaged into a small ellipse whose eccentricity e is given 
by the ratio of the eigenvalues e± of A: 



e 2U/2 



(4.4) 



The vanishing of the Jacobian Det(J) indicates the failure of invertibility of the map ( |4.1| ). 
The loci of points in the lens plane where Det(J) = are called critical lines and the 
corresponding curves in the source plane are called caustics. Critical lines [caustics] 
separate regions corresponding to different numbers of images in the lens [source] plane. 
Therefore, the occurrence of multiple images is signalled by the vanishing o/Det(J), and 
this is the condition that we will study in the following. These concepts are familiar 
from standard gravitational lens theory [^] . We now proceed to study the features that 
are peculiar to lensing gravitational waves. 
Let us consider the Jacobian determinant 

Det(J) = 1 + D - D 2 [d y (5p x ) ■ d x (5py) - d x (5p x ) ■ d y (5 p y)] ; (4.5) 

the divergence d(Sp A ) / dx A vanishes to first order. This follows from the fact that it 
represents the expansion scalar of a congruence of null rays. It was proved correctly 
in Ref. [25|] that this quantity is of second order in the wave amplitudes; we give an 



independent proof in the Appendix. Eq. ( f4.5|) can be written as follows 

Det(J) = 1 + Ji + J 2 , (4.6) 



Ji = ^)D s ^, (4.7) 

J 2 = f(a)D 2 s \d x {5p x ) ■ d y (Spy) - d y (8p x ) ■ d x (5py)\ , (4.8) 

where a = Dls/Ds and f(a) = a 2 {l — a) 2 . The polynomial f(a) is symmetric about 
a = 1/2 (corresponding to Dl = Dls), where it assumes its maximum value 1/16 ~ 
0.0625. 
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From Eq. ( |2.15| ) one obtains the order of magnitude estimate 



J x ~ Dd A (6p A )~^h 2 . 



(4.9) 



Since J2 ~ (D / P) 2 h 2 , we neglect J\ in comparison to J2 in the Jacobian determinant of 
Eq. ( f4.5|) when large values of D/P are considered, which is the case of interest. 

In order for Det(J) to vanish (i.e. to have multiple images of the light source), it 
must be J% + J2 < 0. Since J\ and J2 are terms of order h 2 times D/P or (D/P) 2 , 
respectively (see Eqs. ( |2.12|) and ( fOl) ), in order to have Det(J) = 0, a large value of 
the distance D must balance for the small values of the wave amplitudes. Moreover, for 
large values of D/P, J\ is negligible in comparison with J2, and hence in order to have 
multiple imaging it must be 



f(a) D s d A (5p B ) 



1 



(4.10) 



For standard gravitational lenses, the probability of lensing of a distant source is approx- 
imately maximum when the lens is halfway between the source and the observer p3 



D — 

P 



> 4 - 10 



(4.11) 



This applies also to lensing gravitational waves and therefore values of the numerical 
factor f(a) far from its maximum are not statistically significant. We take f(a) in the 
range - ^ near its maximum. Then, in order to have multiple imaging, the following 
condition must be satisfied: 

5p Dh 
P 

where P is the period of the gravitational wave (computed at the redshift of the lens 
plane, if the D's are chosen to represent angular diameter distances in a FLRW universe) 
and where we used 5p ~ h (see Eq. ( 2.12 )). The inequality ([4.11 ) can be written in the 
form 

(4.12) 



— > S r 



P 



where S c = (4 — VS)c/D. The approximate condition for multiple imaging ( |4.12| ) involves 
the "strength" h and the "size" P of gravitational waves, the geometry of the problem 
(through D) and the fundamental constant c. (|4.12|) is analogous to the well known 



condition for multiple imaging by ordinary gravitational lenses |3T] 



AtiGD 



(4.13) 
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where £ is the two-dimensional (projected or surface) density of the lens and the critical 
density E c depends only from the geometrical factor D and from fundamental constants. 
This is a condition on the "strength" (the mass) and the size of the lens. 

The energy and momentum of high frequency gravitational waves are given by the 
Isaacson effective stress tensor which, in the TT gauge, assumes the form 



T = ° 4 h (TT) h iTT) ) (414) 



32nG 



where ( ) denotes the Brill-Hartle average. For \i — v — 0, T Q0 = pc 2 , where p is the 
effective mass density of the gravitational waves. It is clear from this expression that the 
ratio h/P is a measure of the "mass density" of the particular lens under consideration?]. 
We comment on the analogy between eq. Q4.12| ) and eq. ( [4.131 ): By squaring, dividing 



by G and integrating eq. ( |4.12|) , it is deduced that 

(4.15) 




Strictly speaking, it is eq. ( |4.15| ), not ( |4.12| ) that is the analogous of eq. ( |4.13[ ): the 
quantity £ 9tu has the dimensions of a surface density, and S^Ds/G is the critical surface 
density. 

The particular lens under consideration always produces an odd number of images. 
In fact, the odd image number theorem |35| holds, its proof requiring only boundedness, 
smoothness and transparency of the lens (conditions satisfied by localized gravitational 
waves) . 



5 Comparison between a gravitational wave and an 
ordinary gravitational lens 

It is instructive to compare the action of a gravitational wave with that of an ordinary 
gravitational lens. The latter is a mass distribution described by a Newtonian potential 
$ (satisfying the Poisson equation V 2 $ = 47rp, where p is the lens mass density). It is 
assumed that the lens is smooth, bounded and stationary, i.e 



3 Although the form of the perturbations h^ u depends on the chosen gauge, its order of magnitude, 
and that of h/P, is unaffected by gauge transformations wj). Therefore the approximate condition 



(4.12) for multiple imaging is gauge-independent. 
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• $ is continuous with its first and second derivatives; 

1 /2 

• $ — > and V$ — ► as r = (x 2 + y 2 + z 2 ) 1 — > +oo; 

• d$/dt ~ 0. 

The plane-to-plane map describing the lens action is given by the lens equation and the 
Jacobian matrix can be written |EB] as 



-/i l - x + a 

with 



X = ^, (5-2) 

" s ?/, a S5' (5 - 4) 

where .D and S c have been defined in the previous section and 

/+oo 

dip. (5.5) 
-00 

The Jacobian determinant is given by 

Det(J) = (l- X ) 2 -(A 2 + /i 2 ). (5.6) 

The convergence % describes the action of matter, while A and \i describe the action of 
shear. For a lensing gravitational wave one obtains 

Jgw = ( 1 -m 1 -L ) ' (5 ' 7) 



where 



At = —Dd x (5p x ) = -Dd x (5p x {1) + Spfa) = A gw + 6A 1 , (5.8) 
A 2 = -Dd v (6p y ) = -Dd y (dp y {1) + 5p\ 2) ) = -A gw + 6A 2 , (5.9) 
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fii = -Dd y (5p x ) = -Dd y (5p x {1) + Sp x (2) j = n gw + 5fn , (5.10) 

fi 2 = ~Dd x (5p y ) = —Dd x (dp y {l) + <5pJ 2) ) = fi gw + 5{i 2 • (5.11) 
To first order, one has 

(5.12) 
(5.13) 
(5.14) 
(5.15) 

where the integrals are computed along the photon's path from the source to the observer. 
To first order, A 2 = — Ai = — A gw and /ii = /i 2 = fM gw . In fact, from Eqs. ( |5.12| )-( |5.15| ) 
and from 8a (Sp A ^j = (see the Appendix), it follows that 

D rO , N 



A: = 


I) 

~ ~2 


/ d\ h\\ 
Js 


A 2 = 


D 

~~2 . 


r° 

/ dXh X x' y 
Js 


Mi = 


~ 2~ 


r° 

/ c^A h\\ 
Js 


^2 


D 
~ ~2 


Js 



y ' 



A x + A 2 = -- J s d\ (h xx >\ x + hxx' v , v ) + 0(h 2 

| jf° rfA pV (V% + + 0(/^ 2 ) = + 0(/i 2 ) , (5.16) 



hence A 2 = -Ai + 0(h 2 ). 

In the first order approximation, tensor indices are raised and lowered with 77^" 
and rj^, respectively, and r\ AB = t]ab = 8ab for A, B = 1,2. Hence, to this order, 

h ^' A ,B = h ^' B ,A and A*i = ^2 + 0(/i 2 ). 
The corrections 



= (5p x {2) ) , (5.17) 

5A 2 = -Dd y (6p\ 2) ) , (5.18) 

fyi = -Ddy (5p x {2) ) , (5.19) 

5^2 = -Dd x (6p\ 2) ) , (5.20) 



where Spfys is given by Eq. ( |2.14| ), are of second order in h. To lowest order, we have 

J - = ( 1 -^T i7Z.„)- < 5 ' 21 > 



l 3™ 
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The convergence term is absent and hence the lens action is due only to the shear. This 



result was derived in Ref. [25| using the Raychaudhuri equation and the optical scalars 



formalism, and is now recovered in the vector formalism. 
When the determinant of J gw is considered, one has 

. 2 



Det(J gw ) = 1 - [A 2 gw + /4J + D [d x [Spfa) + d y {5pf 2) )\ . (5.22) 
An order of magnitude estimate as outlined in Sec. 4 gives A 2 , fi 2 ~ (D/P) 2 h 2 >> 



Dd A (5p A ) ~ (D/P) h 2 for large D/P. Therefore, we can write 

Det( J gw ) = 1 - (A 2 gw + fa) (5.23) 

(where the term in brackets is of second order in hD/P). 

It is to be noted that the deflection angle Q2.11D does not depend on the frequency of 
the light. Gravitational waves are achromatic lenses, like ordinary gravitational lenses. 
However, while the latters do not shift the frequency of the photons propagating through 
them, lensing gravitational waves do. The effect is of first order in the wave amplitude 
and has been studied in detail by many authors 0, ||, [7|, || [l|, [14], expecially in conjunc- 
tion with the microwave background anisotropics discovered by the COBE experiment. 
In addition, gravitational waves do not rotate the polarization plane of the electromag- 
netic field, to first order | TDf] . In this aspect, they behave like ordinary gravitational 



lenses. 

It is expected that the images of a distant source created by a gravitational wave 
vary on timescales of the order of the wave period. This could possibly be used to 
explain the variability of some active galactic nuclei or active galaxies. Moreover, the 
details of the images configuration depend on the detailed form and parameters of the 
lensing wave, such as its spatial and time profile, its duration, direction of propagation 
and polarization (see |19[ for an example). 



6 Order of magnitude estimates 

In order to apply the previous theory and the multiple images to be detectable, the 
following conditions must be satified: 

1. geometric optics holds; 
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2. the scale of separation between different images must not be smaller than 10~ 3 
arcseconds. In fact, structures on scales ~ 10 -3 arcseconds can be resolved with 
VLBI, while VLA and optical techniques apply on larger scales^]; 

3. the lens must not be exceptionally rare, i.e. the rate of occurrence of the event 
generating the lensing wave must not be too low; 

4. in order to appreciate variability in the images induced by a lensing gravitational 
wave, its period must not be too short, let us say P < 10 8 s. 

5. to ensure that the point where lensing takes place is in the wave zone of the 
gravitational field, the impact parameter r must obey the condition 

r > \ gw . (6.1) 

In order to satisfy 1), only electromagnetic radiation with wavelength A em satisfying 
X gw ■ (X gw /D) < X em < X gw will be considered. To satisfy 2) note that, if 5 ~ h is the 
deflection angle, it must be 

/i>5-l(T 9 . (6.2) 

3) depends on the particular processes generating gravitational radiation. Since these 
are almost all purely speculative, their rates of occurrence are largely or completely un- 
known and we can only guess their values. A continuous source of gravitational radiation 
will give rise to a permanent lens, while a gravitational wave burst will constitute a tem- 
porary lens. 

The sources of gravitational waves that are most often considered in the literature 

are: 

• stellar collapse with non-spherical symmetry; 

• formation of massive black holes in active galactic nuclei; 

• neutron star collision; 

• black hole collision; 

• close binary systems; 

4 If one is not interested in the possibility of resolving the multiple images created by the gravitational 
wave, but merely in the occurrence of high amplification events, the limit 10 -3 arcseconds can be 
considerably relaxed. 
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• black hole accretion. 



The last two types of objects are continuous sources of gravitational radiation, while the 
others give bursts. In addition, we will consider the stochastic gravitational wave back- 
ground, both primordial or generated ( [|13|, |36| , |37j and references therein). Gravitational 
waves generated by a process involving a body of mass M and size R have dimensionless 
amplitudes (near the source) of order 

hs ~ e — , (6.3) 

where the efficiency e is defined as the fraction of energy radiated away. For processes 
involving neutron stars or black holes one can assume M/R ~ 1/20 and M/R ~ 1, 
respectively. 

Multiple imaging by gravitational waves is, in principle, possible, and it should be 
expected if they satisfy the approximate condition ( f4.12|) . We examine the astrophysical 
sources of gravitational radiation which are most often considered in the literature, in 



conjunction with the condition ( 4.12 ). When the event generating gravitational waves 



involves neutron stars or black holes, the ordinary lensing associated to these objects 



( "microlensing" p3| ) should, in principle, be taken into account: however the separa- 
tion scale between microimages of a distant source created by a compact object is of 
order 10~ 6 arcseconds, not detectable with present techniques, while the multiple im- 
ages due to the gravitational wave may be detectable. In addition, the effects induced 
by gravitational waves vary on a scale different from the typical scale of variation of 
microlensing, and include a frequency shift Su/u ~ h, which is absent in ordinary mi- 
crolensing (and probably unobservable). However, even in the situations in which lensing 
by gravitational waves cannot be separated by ordinary microlensing, the former may be 
dominant, and it is important to study how gravitational waves modify the microlensing 
phenomenon. 

Stellar core collapse 



The research program on collapsing homogeneous ellipsoids by Saenz and Shapiro [[58 



has given the expected maximal efnciences for a "cold" and "hot" equation of state 
as 6 ~ 10 -2 and e ~ 10~ 4 respectively, for a spectrum of emitted gravitational radiation 
broadly peaked between 100 Hz and 1 KHz. Taking the lower value e ~ 10~ 4 one 
gets hs ~ 5 ■ 10 -6 for the wave amplitudes near the collapsing core and the condition 
h ~ hsRs/r > 5 • 10 -9 implies r < 10 3 i?s ~ 10 9 cm. On the other hand it must be 
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f > ^gw — 3 • 10 7 — 10 8 cm; there is a rather narrow permitted range for the impact 
parameter (r ~ 10 8 cm), that gives Dh/P ~ 10 if D ~ 6 • 10 15 cm. 

If the late phase when the ellipsoid has settled down as a rapidly rotating neutron star 
is taken into account, it is found |4]J that the emitted spectrum is very narrowly peaked 
(Au/u ~ KT 3 ) and e ~ 1(T 6 at v ~ 1 KHz, that gives h s ~ 5 • 10" 8 . Condition 
requires r < IO-R5 ~ 10 7 cm; on the other hand, condition (|6.1| ) does not allow for 
lensing on a relevant scale to take place in the wave zone. 

If the core keeps bouncing, its eccentricity becomes large after a sufficient number of 



bounces. This asymmetry |4(| makes the efficiency almost uniformly near its maximum 
value for any initial period above 1 sec to several hundred seconds [}4"2"| |. One gets hs ~ 



1(T 2 at v ~ 1 KHz; the condition (|6~2]) gives r < 2- 10 6 R S ~ 2 - 10 12 cm, while gives 
r > 3 • 10 7 cm. A rather large range of values of the impact parameter is permitted; one 
obtains Dh/P ~ 10 if D ~ 10 12 cm, r ~ 3 • 10 7 cm, or if D ~ 6 • 10 16 cm, r ~ 2 • 10 12 cm. 

Studies of the perturbations of pressureless spherical collapse leading to the for- 
mation of a black hole [||3|-[45] give results that could possibly be extrapolated to 

larger deviations from spherical symmetry | 4"2"|| , getting e ~ 2 ■ 10~ 2 ( j^fs) at v ~ 



IK H z ■ (M /IOMq) 1 for J/M 2 1 (where J is the angular momentum and J = M 2 cor- 
responds to a maximally rotating Kerr black hole). Taking J ~ 0.1 M 2 and M = IOMq, 
it follows that h s ~ 2 • 10~ 6 and (§3) implies r < 400/^ ~ 4 • 10 8 cm, while gives 
r > 3 • 10 7 cm. A rather narrow range of values of r ~ 10 8 is permitted, that gives 
Dh/P ~ 10 if D ~ 3 • 10 16 cm. 

Final decay of a neutron star/neutron star binary 

Rough estimates for the final decay of a binary system composed of two neutron stars 
HI give e ~ 5 ■ 10~ 3 at v < 2 - 3 KHz. Taking v ~ 500 Hz it follows that h s ~ 2.5 • 10" 4 . 
( |6 . 2| ) gives r < 5 ■ 10 4 i?5 ~ 5 • 10 10 cm, while r > 6 ■ 10 7 cm due to ( |6.1| ). A large range 
of values of r is permitted; we get Dh/P ~ 10 if D ~ 3 ■ 10 16 cm, r ~ 5 ■ 10 10 cm, or if 
D ~ 6 • 10 7 cm, r ~ 6 • 10 7 cm. 

Black hole collisions 

The head-on collision of two equally massive, non-rotating black holes has been 
studied numerically ^7j; leading to a single, larger, black hole, with efficiency e ~ 



7-10 4 . If the two initial black holes have nearly enough angular momentum to go into 
orbit before coalescing, a formula derived from extrapolation of perturbation theory 



4"2|, ffffl gives e ~ 3 ■ 10 . This efficiency is expected to hold for P ~ 1 s; (|6.2[) gives 
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r < 6 • 10 6 -Rs ~ 6 • 10 12 cm, while ( |6.1| ) implies r > 3 • 10 10 cm. In the permitted range 
of values of r the condition Dh/P ~ 10 is satisfied if D ~ 3 ■ 10 17 cm, r ~ 3 • 10 10 cm 
and if D ~ 6 • 10 19 cm, r ~ 10 12 cm. 

The case of coalescing black holes was studied also in Ref. [2^]. Although multiple 



imaging was not considered in that paper, the deflection angle of light rays was computed 
using a gauge-invariant formalism. 

The binary pulsar 

The binary pulsar PSR 1913+16 (j|9| and references therein) is believed to radiate 
gravitational waves in a continuous way, according to the predictions of general relativity 
( ||50|| and references therein). The estimated distance of the binary system (believed to 
be a neutron star/neutron star system) is D ~ 5 Kpc and the frequency of the radiation 
is twice the orbital frequency (due to the quadrupole nature of the radiation). From 
these values one obtains^ D/P ~ 3.5 • 10 7 and Dh/P ~ 10 if h ~ 3 • 10~ 7 . An estimate 
of the amplitude of the waves emitted by the binary pulsar gives 

Q MaW 

h ~ — ~ , (6.4) 

r r 

where Q is the quadrupole moment, M is the mass and a ~ 7 ■ 10 10 cm [^TJ] is the 
semimajor axis of the binary system, so that h$ ~ (au/c) 2 ~ 10~ 6 . It is easy to see 



that the conditions (|6.2| ) and (|6.1|) are not compatible. Therefore multiple imaging by 
the gravitational waves emitted by the binary pulsar PSR 1913+16 is impossible. 

The gravitational wave background 



For the gravitational wave background, both primordial or generated (|]r3| and ref- 
erences therein), one has 

P 



h~ ^JVL gw — , (6.5) 

where Po is the present gravitational wave period, R is the radius of the universe and 
Q gw is the cosmological density of gravitational waves (in units of the critical density). 
One has ^ 

-p- ~ \J^Z ^ • (6-6) 



5 The orbital period is ~ 2.8 • 10 4 s M 
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Upper bounds on Q gw have been set in various bands of frequencies. Apart from the 
obvious bound fl gw < 1 for every region of the spectrum, the limit fl gw <C 1 has been 
established in many frequency bands ([TB|, |S(J |57| and references therein). Moreover, 
D/R < 1, hence it is likely that Dh/P <C 1 (the exact value depending on the frequency 
band) and one concludes that multiple imaging by the gravitational wave background 
is, on average, impossible. 



We conclude this section with an example: gravitational waves emitted in a supernova 
collapse in the Virgo cluster. While the dimensionless amplitude hi ~ 10~ 21 on the 
earth (at a distance r± ~ 15 Mpc), one has h,2 ~ h±ri/r2 at a distance T2- If r2 ~ 1 pc, 



h 2 ~ 1.5 • 10~ 4 and assuming P ~ 10~ 3 s, Eq. Q4.10I) tells us that ^ f(a)Dh 2 /P has 
to be greater than, or of the order unity in order to have multiple images. This gives 
(1 — e) _1 < Ds/Dls < where e = 4.4 • 10~ 15 , hence multiple imaging is possible 
with a light source located virtually anywhere beyond the Virgo cluster (of course this 
situation is different if one changes the parameters h 2 , P). 

The amplification on a caustic is infinite in the geometric optics approximation and 
is limited by geometric optics. For ordinary gravitational lenses, A ~ (M/A e m .) 1//3 on a 

1 /2 

fold caustic and A ~ (M/A e . m .) on a cusp caustic, where M is the mass of the lens 



and X e .m. is the electromagnetic wavelength flql . Naively, one expects a similar result 



for lensing gravitational waves, when M is substituted with (h/P) 2 l 3 , where / is the 
size of the gravitational wave packet. However, a realistic computation on the lines of 
p3| is much more complicated and requires the use of precise models that we have not 
considered here. 



7 Discussion and conclusions 



Gravitational waves affect the propagation of light; this is expected for any gravita- 
tional field. The interaction between light and gravitational waves in general relativity 
has been considered by many authors ||— [12|. Recently, the attention was focussed upon 



the contribution of gravitational waves to the microwave background anisotropies dis- 
covered by the COBE experiment [14j], and on the possibility of detecting astrophysically 
generated (as opposed to primordial) gravitational waves through their effects on light 
rays ||19||— 1[22|. The possibility of gravitational waves superimposed to an ordinary grav- 
itational lens was also considered 
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The generalization of the Fermat principle 
to non-stationary spacetimes [0 allows one to approach the problem using a modified 
version of the vector formalism employed for ordinary gravitational lenses p9| . An ap- 



19 



proximate condition (Eq. ( |4.12| )) analogous to that holding for ordinary gravitational 
lenses must be satisfied in order to create multiple images of a distant light source. The 
astrophysical sources of gravitational waves most often considered in the literature have 
been examined. Certain astrophysical sources are shown to satisfy the approximate con- 
dition for the creation of multiple images, due to the balance between the large values 
of the distances involved and the small values of the gravitational wave amplitudes. An- 
other relevant fact is that lensing can take place in regions relatively close to the sources 
of gravitational waves, where the wave amplitudes are larger than those expected in 
the Solar System. The gravitational wave background, instead, is unlikely to produce 
multiple images. As a conclusion, multiple imaging of a distant source by gravitational 
waves is possible in the favourable situations examined, including the collapse of stellar 
cores, the final decay of neutron star/neutron star binaries and black hole collisions. 
Unfortunately, these events (in particular black hole collisions) are very rare and the 
probability of observing the phenomenon is extremely low, mainly because the duration 
of the multiple images is limited to the period of intense emission of gravitational radi- 
ation. Continuous sources of gravitational waves like the binary pulsar PSR 1913+16 
would give a much higher probability, but unfortunately they emit gravitational waves 
too weakly. 

It is to be noted that the limit 5 > 10~ 3 arcseconds in condition 2) of Sec. 6 can 
be considerably relaxed if we do not require that the multiple images created by gravi- 
tational waves be resolved. The high amplification events associated to the presence of 
multiple images would still constitute an interesting and observable phenomenon even if 
the multiple images are not resolved. In fact, this situation occurs in ordinary microlens- 
ing by stars or planets: the typical scale of separation of microimages is 10 -6 arcseconds 
or less, yet the phenomenon has been observed [[51], [52]| . Moreover, it was not necessary 



to consider cosmological distances D in order to have an appreciable probability for 
microlensing: the phenomenon was observed in the Large Magellanic Cloud |5T| and in 
our galaxy |52 |. 

Further work is necessary to assess precisely the probability of observing multiple 
images created by gravitational waves from astrophysical sources, and the details of the 
phenomenon. Contrarily to the conclusions of Ref . [19| , we are not very optimistic on the 
probability of observing such an event, but our results show that this possibility deserves 
some attention. In addition, we restricted our attention to the creation of multiple 
images, while Ref. [l!| include other aspects of lensing. It is also useful to remember 
that, before the discovery of the first gravitational lens system in 1979, gravitational 
lenses were considered mere speculations not occurring in the real world. Today, not 
only the existence of gravitational lenses is universally acknowledged, but gravitational 
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lensing has become one of the most active and promising fields of astronomy. 
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Appendix 

Here we provide an independent proof of the equation 8a(Sp ) = 0, which was found 
also in Ref . . From Eq. (|2.2j ) and from the normalization p^ = we obtain, to first 



order, 

2r)^ 0) 8p v + VP(o)P(o) = > ( A - x ) 

and 

5p° ~^ = \ ( h oo + 2/*o3 + M + 0(h 2 ) . (A.2) 



From Eq. ( pT2l) it follows that d z (5jf) = and Eq. flSJ) gives 

d 3 (/i 00 + 2/i 3 + h 3 ) = + 0(h 2 ) . (A.3) 



Equation ( |2.12| ) gives 5p 3 = and from Eq. ( |A.2| ) it follows that 



Sp° = \ (hoo + 2/io3 + M + 0(h 2 ) . (A.4) 
By setting / = (h 00 + 2/i 03 + h S3 )/2, it follows from Eq. ( |2.12| ) that d 3 f = —d f and 

d 2 00 f ~ dlJ = (A.5) 

along the photon trajectory. Using the linearized Einstein equations for gravitational 
waves propagating outside their sources, 

□ V = > ( A - 6 ) 

d% v = , (A.7) 
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where = h^ u — one obtains 



□ /i 00 = ~Oh 33 = --Uh , Dh 03 = . 



From Eqs. (|A.8|) one obtains 
and 



□/ = 



W) = ofdz (h 0Q + 2h 03 + h 33 )'» = + 0(h 2 
^ Js ,n 

The two-dimensional divergence &a {^p A ^j is given by 
where the last equality follows from Eqs. (|A.5 ) and ( |A.9|) . 



0. 



(A.8) 

(A.9) 
(A.10) 



(A.ll) 
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